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(N 

Q ■ We study generalized Galileons as a framework to develop the most general single-field 

■ inflation models ever, Generalized G-inflation, containing yet further generalization of G- 

' inflation, as well as previous examples such as k-inflation, extended inflation, and new Higgs 

inflation as special cases. We investigate the background and perturbation evolution in this 
model, calculating the most general quadratic actions for tensor and scalar cosmological 
perturbations to give the stability criteria and the power spectra of primordial fluctuations. 
It is pointed out in the Appendix that the Horndeski theory and the generalized Galileons 
are equivalent. In particular, even the non-minimal coupling to the Gauss-Bonnet term is 
included in the generalized Galileons in a non-trivial manner. 

D ■ §1. Introduction 

Scalar fields play important roles in cosmology. On the one hand, inflation in 
the early Universe is now becoming a part of standard cosmology that is driven by 
^ a scalar field called the inflaton.^ 1 ^ The conventional inflaton action consists of a 

£NJ ' canonical kinetic term and a sufficiently flat potential.® [See Ref. 0J) for the latest 

review.] Non-canonical kinetic terms® also arise naturally in some particle physics 
models of inflation such as Dirac-Born-Infeld inflation.® 

On the other hand, it is strongly suggested that the present Universe is dom- 
inated by mysterious dark energy, and its identity might be a dynamical scalar 
field.^ 1 In relation to the present accelerated expansion, modified gravity theories 
have been studied extensively, and in such theories, an extra gravitational degree of 
freedom can often be equivalently described by a scalar field coupled non-minimally 
to gravity or matter. In the decoupling limit of the Dvali-Gabadadze-Porrati brane 
model,® the scalar field has a non-linear derivative self-interaction,® which was later 
generalized to Galileons^® with a number of applications to various contexts in cos- 
mologyP®^23} Thus, in recent years, there have been growing interests in scalar field 
theories beyond the canonical one. 

The most attractive feature of higher derivative theories possessing the Galilean 
invariance d^<p — > d^cf) + is that field equations derived from such a theory contain 
derivatives only up to second orderp® so that it can easily avoid ghosts. Unfortu- 
nately, however, this desired feature ceases to exist once the background spacetime 
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is curved 

To preserve the second-order nature of field equations, the "covariantization" 
of the Galileon has been proposed by Deffayet et alP^ where the theory is no 
longer Galilean invariant. This line of analysis has been further pursued recently to 
yield a more generic class of higher-derivative theories that result in second-order 
field equations.^ In this theory, the Galilean invariance is absent even in the flat 
spacetime limit. 

The purpose of this paper is to provide a comprehensive and thorough study 
of the most general non-canonical and non-minimally coupled single-field inflation 
models yielding second-order field equations making use of Ref . [25]) , which is the most 
general extension of the Galileons but is no longer based on a symmetry argument. 
It would be nice if one could develop a new class of viable inflation models fully 
respecting the Galileon symmetry, (f> — >■ <f> + b^x^ 1 + c, as was attempted in Refs. [2Tj). 
[22]) , but such a symmetry must be actually broken to construct a phenomenologically 
viable inflation model, namely, to terminate inflation and reheat the Universe. Thus, 
our strategy here is more similar to the G-inflation modefS} and is indeed the most 
general extension of it. Hence, one may call the model presented here as Generalized 
G-inflation or G 2 -inflation. 

Special cases of the generalized Galileons can be derived from a relativistic probe 
brane embedded in a five-dimensional bulkj^"^ and hence, they are possibly re- 
lated to fundamental theory and particle physics. The scalar field theories we are go- 
ing to study thus include not only all the previous examples considered in the context 
of single-field inflation, but also recent developments and their further generaliza- 
tion. We clarify the generic behavior of the inflationary background and investigate 
the nature of primordial tensor and scalar perturbations at linear order. Given a 
specific model, our formulas are helpful to determine the evolution of cosmological 
perturbations and its observational consequences. 

This paper is organized as follows. In the next section, we define the scalar-field 
theories that we consider. We then provide the background cosmological equations 
in §3, and explore two possible inflationary mechanisms. In §4, cosmological pertur- 
bations are considered and the quadratic actions for tensor and scalar perturbations 
are computed, which are used to give stability conditions and evaluate the primordial 
power spectra. Our conclusion is drawn in §5. 

§2. Generalized Higher-order Galileons and Kinetic Gravity Braiding 

Galileons^ and their covariant extension^ have been further generalized re- 
cently to yield the most general scalar field theories having second-order field equa- 
tions.^ 1 The first two terms of the generalized Lagrangian corresponding to (deft) 2 
and (dcp^Ocj) in the original theory are given by^'^ 



C 2 = K(<j>,X), 

C 3 = -G 3 ((j>,X)n^ 



(2-1) 
(2-2) 
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where K and G3 are generic functions of 4> and X := —d^cj)d^(f)/2. Similarly, higher- 
order Galileons can be generalized to give^J 

£ 4 = G 4 (0, X)R + G AX [(n^) 2 - (V^V^) 2 ] , (2-3) 

£5 = G 5 (0, XJG^VV^ - ^f- [(D0) 3 - 3 (D0) (V^V^) 2 + 2 (V^V^) 3 ] ,(2-4) 

where is the Ricci tensor, G^ is the Einstein tensor, (V^V^) 2 = V^V^V^V^, 
(V M V„</>) 3 = V^V^VVVVaV^, and G iX = dG l /8X. Setting G 3 = X, G 4 = 
X 2 , and G5 = X 2 , the above Lagrangians reproduce the covariant Galileons intro- 
duced in Ref . |24"|) . The non- minimal couplings to gravity in £4 and £5 are necessary 
to eliminate higher derivatives that would otherwise appear in the field equations. 
Note that we do not need a separate gravitational Lagrangian other than £4; for 
G4 = Mpj/2, £4 reduces to the Einstein-Hilbert term. We also obtain a non-minimal 
coupling of the form f((f>)R from £4 by taking G4 = f((j>)- The non-standard kinetic 
term G Ml 'd^cj)d u 4> that is considered, such as in Ref. I30|) . turns out to be a special case 
G5 oc (j) of £5 after integration by parts. Equation (24) of Ref. I3T|) . which is obtained 
from a Kaluza-Klein compactification of higher-dimensional Lovelock gravity, turns 
out to be equivalent to £5 with G5 = — 3X/2. 

We thus consider a gravity + scalar system described by the action 



5 

/ d 4 xV=gjCi, (2-5) 



S = 

i=2 

which is the most general single scalar theory resulting in equations of motion con- 
taining derivatives up to second order. This action contains only four independent 
arbitrary functions of <f> and X. This theory represents a general class of single-field 
inflation, including models that have not been studied so far, as well as almost all the 
previously known models such as potential-driven slow-roll inflation,® k-inflationP 
extended inflation,^} and even new Higgs inflation® a s special cases 1*^1 
Note in passing that 

XR + (Ocj)) 2 - (V^Vufj)) 2 = G^V^fiVcp + total derivative 

= -^.VV^ + total derivative, 

which implies that £4 with G 4 = X and £5 with G5 = — <f> are equivalent. Similarly, 
£3 with G3 = f(4>) and £2 with K = —2Xf ( j > coincide up to total derivative. These 
facts can be used to check the calculations. 



§3. Background equations 



Let us derive the equations of motion describing the background evolution from 
(|2 5j) . The easiest way is to substitute <p = <f>(t) and the metric ds 2 = —N 2 (t)dt 2 + 

*' Even the curvature-square inflation model as well as more general f(R) inflationP a which 
do not contain any scalar field and result in fourth-order equations of motion, can be recast in the 
present form by defining a new field as (j> = d//d_R. 
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a 2 (i)dx 2 to the action. Variation with respect to N(t) gives the constraint equation 
corresponding to the Friedmann equation, which can be written as 

5 

£> = 0, (3-1) 

where 

£ 2 = 2XK X - K, (3-2) 

£ 3 = 6X<PHG 3X - 2XG 3(f> , (3-3) 

£ 4 = -6H 2 Gi + 24H 2 X(G 4X + XG 4XX ) - V2HX^G HX - GH^G 4(f> , (3-4) 

S 5 = 2H 3 X<j> (5G 5X + 2XG 5XX ) - 6H 2 X (3G b<j) + 2XG b4>x ) . (3-5) 

The above quantities contain derivatives of the metric and the scalar field up to first 
order. 

Variation with respect to a(t) yields the evolution equation, 

5 
i=2 

where 

V 2 = K, (3-7) 

V 3 = -2X (G 3(j> + 4>G 3X ) , (3-8) 



V A = 2 ( 3H 2 + 2H ) G A - \2H 2 XG AX - 4HXG 4X - 8HXG 4X - 8HXXG AXX 



+2 (0 + 2H4>) G 4<p + ULG m + 4X (0 - 2H<i>) G 4<px , (3-9) 
V 5 = -2X (2H 3 + 2HH0 + 3H 2 ^j G 5X - 4H 2 X 2 4>G 5XX 

+4HX (X - HX) G 5<px + 2 [2 (HX)'+ 3H 2 X] G b<f> + AHX^G^. (3-10) 



The background quantities £{ and V% are defined in an analogous way in which 
the energy density and the isotropic pressure of a usual scalar field are defined. In 
the present case, however, the distinction between the gravitational and scalar-field 
portions of the Lagrangian is ambiguous, and hence, in that sense, the gravitational 
contribution is included in the above expressions. Indeed, one can see that, for G A = 
Mpj/2, £ 4 and V A reduce to the "minus" of the Einstein tensor G£\ £4 = -3M^H 2 
andP 4 = M 2 X {3H 2 + 2H). 

Variation with respect to 4>(t) gives the scalar-field equation of motion, 

where 

J = 4>K X + 6HXG 3X - 24>G 34> + m 2 4> {G AX + 2XG AXX ) - l2HXG A<t>x 

+2H 3 X (3G 5X + 2XG 5XX ) - 6H 2 <j> (G 50 + XG^ X ) , (3-12) 
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and 

r 2 



P^ = K 4> - 2X (Gw + <j>G^x) + 6 [2H Z + H) G H + 6H (X + 2HX) G H x 

-6H 2 XG 5H + 2H 3 X<j>G 5<px . (3-13) 

Not all of the equations (|3-1|) . (|3 6|) . and (|3-lip are mutually independent: Eq. (|3-11|) 
[ d£SD ] can be derived from Eqs. (03} [ (pMB ] and (l34|) . 

Below, we present two extreme cases of inflation, one driven purely kinetically, 
which is an extension of the original G-innatiorPJ corresponding to the case with 
G4 = Mpj/2 and G5 = 0, the other driven by a scalar potential as an extension of 
Higgs G-inflationP3 

3.1. Kinetically driven G-inflation 

Let us start with a shift-symmetric model, (ft — > (ft + c. This in particular implies 
that 4> does not have any potential. In this case, the field equations are 

5 

£i = <P J - K - 6H 2 (G 4 - 2XG AX ) + 4H 3 Xcj)G 5 x = 0, (3-14) 

i=2 

5 , 

J2(£i + Vi) = 4>J- 2X4>G 3X + 2- [2H(G A - 2XG 4X ) - H 2 X<j)G 5X 

i=2 

= 0, (3-15) 
^(a 3 J) = 0. (3-16) 

From Eq. (|3T6p . one immediately finds that J oc a~ 3 — > 0. Shift-symmetric models 
thus have an attractor, J = 0, along which H =const, =const, satisfying 

<PK X + 6HXG 3X + 6H 2 <j) {G iX + 2XG iXX ) 

+2H 3 X (3G 5X + 2XG 5XX ) = 0, (3-17) 
K + 6H 2 (Gi - 2XG 4X ) - 4H 3 X(j)G 5X = 0. (3-18) 

Provided that Eqs. (|3- 17|) and (|3- 18|) have a non-trivial root, H 7^ 0, 4> 7^ 0, we 
obtain inflation driven by (fi's kinetic energy. This is the generalization of kinetically 
driven G-inflatiorP^J 1 and kinetic gravity braiding.^ 

Although the shift-symmetric Lagrangian can nicely accommodate a de Sit- 
ter solution as an attractor, the shift symmetry must be broken in some region in 
the field space to end inflation and to reheat the Universe. Following the argu- 
ments in Refs.P'ES a graceful exit from kinetically driven inflation is possible with 
gravitational reheating.® However, a detailed analysis of the reheating stage after 
kinetically driven inflation is beyond the scope of the present paper. 

3.2. Potential- driven slow-roll G-inflation 

Suppose that the functions in the Lagrangian can be expanded in terms of X as 

K(4>,X) = -V{(j>)+lC{(j>)X + -- - , (3-19) 
Gi(4>,X) = 9i {4>) + hi(<t>)X + • • • , (3-20) 
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and consider the case in which the inflaton field value 4>(t) changes very slowly. In 
this case, the potential term manifestly breaks the shift symmetry, and thereby, 
the model is capable of a graceful exit from inflation.® Neglecting all the terms 
multiplied by 4> in the gravitational field equations, we obtain 

5 5 

Y, V i--Y. £i -- V ^ + Q 3^)H\ (3-21) 

i=2 i=2 

where we have assumed 

\H\ < H 2 and \4>\ < \Hj>\. (3-22) 
We may thus have slow-roll inflation with 

H 2 ~ (3-23) 

During slow-roll, we approximate 

\j\ < \HJ\, \gt\ < \H 9i \, \hi\ < \Hhi\. (3-24) 
Under the above approximation, we have the slow-roll equation of motion for <f>, 

3HJ~ -V^ + l2H 2 g Hl (3-25) 

with 

J ~ K<j> - 2g 3(f> <p + 6 (HhaX + H 2 h^ - H 2 g b4> ^ + H 3 h 5 X^j . (3-26) 

Which term is dominant in Eq. (|3-26p depends on the magnitude of the coefficients 
hi((f)) of X. Note here that we can set (73 = and 95 = without loss of generality, 
because g^ can be absorbed into the redefinition of K. and g^ into /14, that is, 
£ - 2d3cj> -*1C, h^ - g^ — > hi. 

Equations (|3-23|) and (|3-26|) imply that the dominant contribution to the in- 
flationary Hubble parameter is the potential V in £2, while any of the terms in 
Eq. (|3-26[) can participate to determine the actual dynamics of the scalar field. 
Therefore, the slow-roll parameters expressed in terms of the potential may look 
very different from the standard ones in general. This is the generalization of the 
Higgs G-inflatiorPJ (see also Ref.ESD). 

§4. Quadratic actions for tensor and scalar perturbations 

In this section, our goal is to compute quadratic actions for tensor and scalar 
cosmological perturbations in Generalized G-inflation. We use the unitary gauge in 
which (j) = <f)(t) and begin with writing the perturbed metric as 



ds 2 = -N 2 dt 2 + 7y (dx i + N'dt) (dx j + N j dt) , 



(4-1) 
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where 



N 



1 + a, N i = ftp, Jij = a 2 (t)e 2C + hij + ^h ik h k ^j . (4-2) 

Here, a, /3, and £ are scalar perturbations and is a tensor perturbation satisfying 
ha = = /iijj. With the above definition of the perturbed metric, \J—g does not 
contain up to second order, and the coefficients of C, 2 and vanish, thanks to 
the background equations. 

4.1. Tensor perturbations 

The quadratic action for the tensor perturbations is found to be 



S, 



(2) 
T 



1 



dtd 3 xa 3 



where 

F T := 2 [G A - X ^G 5X + G 5(j> 
Qt := 2 G4 — 2XG±x — X (^HcfyG^x — G^ 
One may notice that Qt can also be expressed as 



Qt 



1=2 



The squared sound speed is given by 



2 _ F T 

C T — 

SdT 



(4-3) 

(4-4) 
(4-5) 

(4-6) 
(4-7) 



One sees from the action (|4-3p that ghost and gradient instabilities are avoided 
provided that 



F T > 0, Q T > 0. 



(4-8) 



Note that c\ is not necessarily unity in general cases, contrary to the standard 
inflation models. 

To canonically normalize the tensor perturbation, we define 



d^T := —dt, z T := ^ (F t Qt) 
a 2 



1/4 



ZThij, 



and then the quadratic action is written as 



S. 



(2) 



1 



T - 2 I dVTd X 



J \2 



(Vvij) 2 + -^v 



T„,2 



ZT 



(4-9) 



(4-10) 



where a prime denotes differentiation with respect to i/t- In terms of the Fourier 
wavenumber k, sound horizon crossing occurs when k 2 = z'J^/zt ~ V^r f° r each 
mode. 
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On superhorizon scales, the two independent solutions to the perturbation equa- 
tion that follows from the action (j4- lOj) are 

oc zt and zt f — ?g-- (4 - H) 

In terms of the original variables, the two independent solutions on superhorizon 
scales are given by 

hij = const and / . (4 , 12) 

J a 6 g T 

The second solution corresponds to a decaying mode. 

To evaluate the primordial power spectrum, let us assume that e := —H/H 2 ~ 
const, 

/t := — const and gr '■= ~ const. (4T3) 

We also define the variation parameter of the sound velocity of tensor perturbations 

as 

ST '-=W^=2 {fT ~ 9T) - (4 ' 14) 

Clearly, only two of the three parameters are independent. We additionally impose 
conditions 

1 - e - / T /2 + g T /2 > 0, (4-15) 

3-e + g T >0. (4-16) 

The former (equivalent to 1 > e + st) guarantees that the time coordinate y? runs 
from — oo to as the Universe expands. The latter implies that the second solution 
in (|4-12p indeed decays. We see that zt can be written as 

ZT = ^L 1 [(-^)/(-^)i 1/2 - (4 . 17) 

where the quantities with * are those evaluated at some reference time yx = yr*- 
The normalized mode solution to the perturbation equation is given in terms of the 
Hankel function: 

rW 



yH^{-ky T )e ih (4-18) 



7T 



2 

where 



VT := lzl±g , (4 . 19 ) 
2 - 2e - } T + gr 
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and eij is a polarization tensor. Notice that the conditions (|4T5p and (|4-16p guar- 
antee the positivity of vt- On superhorizon scales, — kyx C 1, we obtain 



2 VT 



- 2 r{ur) (-yr) 1/2 - UT ,j 3/2 - VTi 



r(3/2) z T 

Thus, we find the power spectrum of the primordial tensor perturbation: 



Vt = 87T 



G 1 ' 2 H 2 



x-3/2 4^-2 



(4-20) 



(4-21) 



-ky T =l 

where 7 T = 2 2 ^~ 3 |r(> T )/T(3/2)| 2 (l - e - / T /2 + gr/2). The tensor spectral tilt is 
given by 



tlt = 3 — 2z^. 



(4-22) 



Contrary to the predictions of the conventional inflation models, the blue spec- 
trum riT > can be obtained if the following condition is satisfied, 



4e + 3/ T - 9t < 0. 



(4-23) 



This condition is easily compatible with the conditions (|4- 15|) and (|4-16p . Thus, 
positive and large qt compared with e and Jt can lead to a blue spectrum of tensor 
perturbations. In deriving the above formulas, we only assumed that e, fx, and qt 
are constant. These parameters may not necessarily be very small as long as the 
inequalities (|4-15p and ()4-16p are satisfied under a sensible background solution. 

4.2. Scalar perturbations 

We now focus on scalar fluctuations putting hij = 0. Plugging the perturbed 
metric into the action and expanding it to second order, we obtain 



S 



(2) 



dtd 3 xa d 



v 2 -v 2 • v 2 ■ 

-20a— P + 2QtC^P + 60«C - 2G T a-^C 



a- 



a- 



(4-24) 



where 



E := XK X + 2X 2 K XX + \2R4>XG 3X 
+6H<j>X 2 G 3 xx - 2XG M - 2X 2 G 



6H G4 



+6 



H 2 (7XG AX + lQX 2 G 4X x + 4X 3 G AX xx) 



—H<j) (G40 + fiXG^x + 2X 2 G^xx ] 
+30H 3 <j)XG 5X + 26H 3 <j)X 2 G 5X x 



(4-25) 



:-- 



+4H 3 <pX 3 G 5X xx - 6H 2 X(6G 5(P + 9XG 5 ^x + 2X 2 G^ X x), 
-<pXG 3X + 2HG 4 - mXG iX - 8HX 2 G± XX + <j>G 4(P + 2X^G 44X 
-H 2 ^ {5XG 5X + 2X 2 G 5X x) + 2HX (3G 5 ^ + 2XG b<f>x ) . (4-26) 
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It is interesting to see that even in the most generic case, some of the coefficients 
are given by Tt and Qt, i.e., the functions characterizing the tensor perturbation, 
and only two new functions show up in the scalar quadratic action. Note that the 
following relations hold: 



E 







i=2 
1 5 

-Y 



i=2 



dH 



6 



i=2 



d£i_ 



(4-27) 
(4-28) 



which compactify the above lengthy expressions. 

Varying the action (|4-24p with respect to a and /3, we obtain the constraint 
equations 



V 2 • V 2 



0. 



9a - g T c = 0. 



(4-29) 
(4-30) 



Using the constraint equations, we eliminate a and (3 from the action (|4-24|) and 
finally arrive at 



S 



(2) 



GsC 2 



where 



?s ■= ~ ( ^Qt 



1 d / a 
adl \0 { 

Qs ■= ~q2~Gt + 3Qt- 



Ft, 



(4-31) 

(4-32) 
(4-33) 



The analysis of the curvature perturbation hereafter is completely parallel to 
that of the tensor perturbation. The squared sound speed is given by c| = J-s/Qs, 
and ghost and gradient instabilities are avoided as long as 



F s > g s > 0. 



(4-34) 



In the case of k-inflation where G3 = = G5 and G4 = Mpj/2, we have J~s = 
Mpje. This implies that the interesting regime H > is prohibited by the stability 
requirement in k-inflationP3 However, the sign of H and the stability criteria are 
not correlated in more general situations. This point was already clear in G-inflation 
and kinetic gravity braiding for which G3 ^ oP^'^S Stable cosmological solutions 
with H > offer a radical and very interesting scenario of the earliest Universe ED'EHl 1 

The stability conditions (|4-34p for the scalar perturbation as well as (|4-8|) for the 
tensor perturbation have been derived in the case of the covariant Galileon, for which 
K = -c 2 X,G 3 = -c 3 X/M 3 ,Gi = M^/2 - c 4 X 2 /M 6 , and G 5 = 3c 5 X 2 /M 9 PJ It 
can be verified that our general formulas correctly reproduce the result ofP^f 
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Using the new variables 

dy s :=^dt, z s :=V2a(F s Qs) ll \ u := Z S (, (4-35) 
a 

the curvature perturbation is canonically normalized and the action is now given by 



4 2) = \ I d Vsd 3 x 



( u ')2 _ (y u f + z lj 
zs 



(4-36) 



where a prime denotes differentiation with respect to ys- Each perturbation mode 
exits the sound horizon when k 2 = z'L/ ' z$ ~ where k is the Fourier wavenumber. 

The two independent solutions on superhorizon scales are 

C = const and / » . (4-37) 
J a 6 Gs 

During inflation, it may be assumed that Q$ is slowly varying. In this case, the 
second solution decays rapidly Note, however, that the non-trivial dynamics of the 
scalar field can induce a temporal rapid evolution of Qg, which would affect the 
superhorizon behavior of the curvature perturbation through the contamination of 
the second mode in the same way as in Ref. [39]) . Given the specific background 
dynamics, one can evaluate such an effect using our general formulas. 

Closely following the procedure we did in the case of the tensor perturbation, 
we now evaluate the power spectrum of the primordial curvature perturbation. To 
do so, we assume that e ~ const0 

fs ■= -TTE~ - const ' 9S ■= -fttt - const > ( 4 ' 38 ) 



and then define 



Us := o — o — ' 4-39 

2 - 2e - fs + gs 



The power spectrum is given by 



C 2 3/2 47r 2 



(4-40) 



-ky s =l 

where 75 = 2 2!ys ~ 3 |r(i/ 5 )/r(3/2)| 2 (l - e - f s /2 + g s /2). The spectral index is 

n s - 1 = 3 - 2u s . (4-41) 
An exactly scale-invariant spectrum is obtained if 

e + \fs ~ \gs = 0. (4-42) 



*^ By defining the variation parameter of the sound velocity of scalar perturbations as ss 
cs/(Hcs) = (fs — gs) /2, the formulae with fs and/or gs can be rewritten in terms of cs- 
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Here again, e, fs, and gs are not necessarily very small (as long as n s — 1 ~ 0). 
Taking now the limit e, /t,9t, fSi9S *C 1, the tensor-to-scalar ratio is given by 

Note that even in the de Sitter limit where e, Jt,9t, fs,9S 0, the scalar pertur- 
bation can be produced in general, r / 0. 

In the case of potential-driven slow-roll inflation in §3.2^ we have, to leading 
order in slow-roll, 

?S*-^{!C + GH 2 h 4 ) +^(h 3 + H 2 h 5 ) , (4-44) 

& - §2 + QH2h *) + fa + , (4-45) 

and J~t — Qt — 2g±, where we used the slow-roll equation 2g±e + g^/H ~ /2H 2 . 
In this case, we have ~ 1 and n-r ~ — (2e+g-r) with fx~gr — 9i/{Hg±). If the /C 
or /14 term dominates in J, we have c| ~ 1 and T s - Gs - J<j>/(2H 2 ) ~ 54 (2e + g T ), 
which yields the standard consistency relation: 

r ~ -8n T . (4-46) 

On the other hand, if the or hs term dominates, then we have c s — 2/3, J-g ~ 
2J<p/{3H 2 ) ~ (4 5 4/3)(2e + 5T), and - J<p/H 2 ~ 2# 4 (2e + c/ T ), which yields a new 
consistency relation!^ 1 

r ~ Q~ nT ' (4-47) 

Thus, we can discriminate which term dominates in the dynamics using the consis- 
tency relations. 



§5. Summary 



In this paper, generic inflation models driven by a single scalar field have been 
studied. Our gravity + scalar-field system is described by the generalized Galileons, 
which do not give rise to higher derivatives in the field equations despite the non- 
minimal coupling, e.g., of the form G{<p,X)R. The class of inflation models is the 
most general ever proposed in the context of single-field inflation. 

We have seen that if the Lagrangian has a shift symmetry, <j> — > <fi + c, de Sitter 
attractors are present and, hence, inflation can be driven by $>'s kinetic energy. 
Reheating after kinetically driven inflation is possible by breaking the shift symmetry, 
but the way to break it depends on the explicit construction of the originally shift- 
symmetric Lagrangian itself. 

We have also derived slow-roll equations of motion for potential-driven inflation, 
in which the scalar-field dynamics is modified by the higher-order Galileon terms. 

We have determined the most generic quadratic actions for tensor and scalar 
cosmological perturbations. Using them, we have presented the stability criteria for 
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both types of perturbations. The primordial power spectra have also been computed. 
Note that, since the propagation speeds of the two types of fluctuations can be 
different, we must evaluate the power spectra for the same comoving wavenumber 
at different epochs, which may have some consequence 

It would be interesting to extend the present linear analysis of the curvature per- 
turbation to non- linear order along the line of Refs. l23|) J4ip . In relation to cosmolog- 
ical perturbations beyond linear order, it would be important to evaluate primordial 
non-Gaussianities generated from Generalized G-inflation. 

Acknowledgements This work was supported in part by JSPS Grant-in-Aid 
for Research Activity Start-up No. 22840011 (T.K.), Grant-in-Aid for Scientific 
Research Nos. 23340058 (J.Y.) and 21740187 (M.Y.), and Grant-in-Aid for Scientific 
Research on Innovative Areas No. 21111006 (J.Y.). 

Appendix A 

The Horndeski action, generalized Galileons, and non-minimal coupling to 

the Gauss-Bonnet term 

In 1974, Horndeski presented the most general action (in four dimensions) con- 
structed from the metric g^, the scalar field <fi, and their derivatives dg^ u , d 2 g liu , d^g^, 
■■■ ,d(ft,d 2 (j),d S (ft,---^b still having second-order field equations. The Horndeski 
theory has been revisited recently by the authors of Ref. I4"3|) . In this appendix, we 
point out that the Horndeski theory and the generalized Galileons are equivalent. 

In terms of the notation of Ref. |4U|) but using X = —d^d^fyjl rather than 
p = d^(pd^(j), the Lagrangian of the Horndeski theory is given by 



£ §afa 



+2K 3 xV a </>V / VV l/ V / #V <T V, 



+ 8$ [(F + 2W)RJ V + 2F x V^V a cpV u V^ 



+2K S V a (jN ll (IN v V^\ — 6 (F0 + 2W^ - Xk 8 ) + Kg, (A-l) 

where 5^^ 2 2 ';;^™ = n\8^8?£...8^} , and k±, k 3 , k 8 , and Kg are arbitrary functions of 
<j) and X. We also have two functions F = F((p,X) and W = W(4>), and the former 
is constrained so that Fx = 2(k 3 + 2Xk 3 x — ^i</>), while the latter can be absorbed 
into a redefinition of the former. We are therefore left with four arbitrary functions 
of 4> and X, in accordance with the generalized Galileon. 

The above Lagrangian can be mapped to that of the generalized Galileon by 
identifying 

K = K 9 + iX / AX' {k^ - 2k 3H ) , (A-2) 

rX 

G 3 = - 2Xk 8 - 8Xk 3 ^ + 2 / dX'(K 8 - 2k 3 ^), (A-3) 

G 4 = 2F- 4Xk 3 , (A-4) 
G 5 = -4asi, (A-5) 
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where we redefined F so that F + 2W — > F. Now we see that the two theories are 
in fact equivalent. In deriving the Lagrangian (|A-ip . Horndeski started from the 
assumptions that are weaker than those made by Deffayet et al.pS) although the 
latter worked in arbitrary dimensions. 

Since the generalized Galileon is the most general theory in four dimensions com- 
posed of g^ u , 4>, and their derivatives, which gives the second-order field equations, 
it must reproduce the non-minimal coupling to the Gauss-Bonnet term,® 

C(cP) {R 2 - AR„ V R^ + R^R^p°) , (A-6) 

which seems non-trivial at first glance. One can show that, by taking 

K = 8£ (4) X 2 (3-lnA), (A-7) 

G 3 =4£ (3) X(7-31nX), (A-8) 

G 4 = 4^X (2-lnX) , (A-9) 

G 5 = -4£ (1) lnX, (A- 10) 

or, equivalently, k\ = ti^'lnX, K3 = £( 2 )lnX, Kg = 0, and Kg = 1G£^X 2 , where 
^(n) ._ Qn^jQfj/h^ generalized Galileon indeed reproduces the non- minimal cou- 
pling of the form (|A-6j) . Probably the shortest way to confirm this fact is to substi- 
tute Ki to the field equations presented in Ref. [4~2|) and to compare them with those 
obtained from (|A-6j) . 

Similarly to f(R) gravity, the gravitational theory described by 



£ = f + f(&), &--=R 2 - W^RT + R^paBT^, 



where f(&) is an arbitrary function of the Gauss-Bonnet term, contains an extra 
scalar degree of freedom, and hence, (|A-11|) must be recast in the Lagrangian of the 
generalized Galileon. Noting that the Lagrangian (|A-llj) can be equivalently written 
as 

£ = f + f(<i>) + f <!>(&- 4), (A-12) 

and the non-minimal coupling fJS is reproduced by Eqs. (|A-7p ~ (|A-10p . it is now 
straightforward to translate (jA-lip to the generalized Galileon. 

It is easy to see explicitly in the cosmological equations of motion that the 
contribution from the non-minimal coupling (|A-6P can indeed be reproduced from the 
non-trivial functions (|A-7p - (|A-10p . Both the generalized Galileon with (|A-7j> — (|A- 10|) 
and the Lagrangian (|A-6|) give the following identical contributions to the background 
and perturbation equations: for the background gravitational field equations, 

£ D -2AH 3 i, (A- 13) 



H 2 '£ + 2H [H 2 + H)£ 



(A-14) 



V D 

for the background equation of motion for <fr, 

D 2AH 2 ( H + H 2 ) fy, (A-15) 
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and for the quadratic actions of the tensor and scalar perturbations, 
T r D 8& Qt D 8Hi E D -48# 3 £, O d 12^^. 



(A-16) 



Appendix B 

FieW equations 



In this Appendix, we present both gravitational- and scalar-field equations de- 
rived from the action (|2-5|) for completeness Varying the action, we obtain 



where 



7 [LV 



7 [LV 



-g 



i=2 



E + E W - yMJ , 

-(-total derivative 



i=2 



1 



1 



GaG^ u - ^GixRV^Vrf - \g axx [(Q0) 2 - (V^V^) 2 ] V^V^ 

-G 4 xn</»V M V^ + G 4X V A V^V A V^ + 2V A G 4X V A V (A >V i , ) </> 
-V A G 4 jV A ^V^ + ^ (G^D0 - 2XG iH ) 



+9nu< -2G 4(/)X V Q V^^V a ^V /3 ^ + G4X* VaV A 0V0V A <£V a 0V 



+-G 4X [(°<A) 2 - (V c * 



+ 2 



9 (XV 



G 4X i? Q/3 V 



V A G 4X V A (/>D< 



-V (/1 G 4 xV l/) ^D^ 

+G^ x R IU3a , p V a <l>V p <f> - G^V^V^ - G^V^V^ 
+2G 4(/ , X V A 0V A V ( ^V^ - G 4X xV a (/)V a V At 0V^0V^V l 



7 X rh\7°"Ul 3 A ' ^ -D V7 V7A, 



~2 v 0* [GsxV ^] V a V v) <l>n<l> + -V (jU [Gg^V^] 
-V A [G^V^] V^vV 
+^ [v A (g 50 V a </>) - V Q (G 5X 



v a v 



V (M G 5 G ! , )A V A ( 



+V a G 5 V^ J R a( ^ )/3 
+ ^V (/i G 5 xV^ [(a</>) 2 - (V ' 



(B-l) 

(B-2) 
(B-3) 



(B-4) 



*' In the course of the preparation of this manuscript, we became aware that Gao has also 
calculated gravitational field equations in the present modelP^ After some iterations, his result has 
converged with ours and we are in full agreement. 



16 



T. Kobayashi, M. Yamaguchi, and J. Yokoyama 



1. 



V, 



+-V>V a G 5 x 
-^GsxG^V^V^V^ - ^G 5X ncpv a V^V a V^ 
+lG 5 x(n^) 2 V M V^ + ^G 5 xx [(n^) 3 - 3Q0(V Q V^ 
+2(V Q V / #) 3 ] V^V^ + ^V A G5G MI/ VV 



SD^V^V^) 2 + 2(V Q V^) 3 ] + V a G 5 # 



1. 



-^V a (G 50 V» □<£ + (G 50 V^) V^V^ - -V Q G 5X V Q ID( 
+^V Q G 5 iV^V Q V^ - iv A G 5 xV A [(□(A) 2 - (V a V^) 2 ] 



(B-5) 





^ = fy, 




(B-6) 




Pi = V„G 34> V»<i>, 




(B-7) 




P% = G 44 ,R + G HX [(□'A) 2 - (V M V^) 2 ] , 




(B-8) 




p| = -v^crvu^ - \g^ x [(^f - 


7 M V^) 2 + 2(V M V^ 


O 3 ] ,(B-9) 


and 








J 2 - 






(B-10) 


J 3 - 


-C 3X y^ + G 3 xV M X + 2G 3 ^V^, 




(B-ll) 


J 4 - 


->C4xV M + 2G 4 xiV v > " 2G 4xx (n^X - 








-2G 4( px(0(f>V^ + V^X), 




(B-12) 


J 5 - 


-j0 5X V^ - 2G 50 G^V l > 







-G5X 



+G 5X x {^V M X [(a</>) 2 - (V^V^) 2 ] - (□<AV M V> - V a V^V Q V»j 
+G^x j^V^ [(D0) 2 - (VaV^^) 2 ] + D0V M X - VXV^V^j . (B-13) 



The gravitational- and scalar-field equations are thus given by 

XX, = 0, ^(j2Ji)=i:n- 



(B-14) 



i=2 



,i=2 



j=2 



One might worry that V M gives rise to higher derivatives as apparently contains 
second-order derivatives. However, this is not the case because the commutations of 
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higher derivatives can be replaced by the curvature tensors and thus are canceled. 
For instance, one has V M (D^W + VI) = (D^) 2 - (V a Vp<j)) 2 - R^V^cpVcj). 
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